Selective population and neutron decay of the first excited state of semi-magic 23 O 
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We have observed an excited state in the neutron-rich semi-magic nucleus 23 O for the first time. 
No such states have been found in previous searches using 7-ray spectroscopy. The observation 
of a resonance in n-fragment coincidence measurements confirms the speculation in the literature 
that the lowest excited state is neutron unbound and establishes positive evidence for a 2.8(1) MeV 
excitation energy of the first excited state in 23 0. The non-observation of a predicted second excited 
state is explained assuming selectivity of inner-shell knockout reactions. 
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Magic numbers are pillars of our understanding of nu- 
clear structure. For each magic number of nucleons one 
can observe, e.g., (i) a large number of stable isotopes or 
isotones, large binding energies, and corresponding struc- 
tures in separation-energy systematics, (ii) increased ex- 
citation energies of first excited states for even-even nu- 
clei, and a corresponding reduction of their ground-state 
transition matrix elements, and (iii) a decreased level 
density at low and moderate excitation energies 0]. Nu- 
clei with one nucleon added or removed from a doubly- 
magic core are well described by a single-particle (SP) 
model; their properties such as spins and parities, mag- 
netic moments, energies of the lowest excited states, and 
transition matrix elements are well reproduced by SP es- 
timates. The underlying physics for the emergence of 
magic numbers is the presence of gaps in the SP spec- 
trum around the Fermi energy, which makes the last nu- 
cleon tightly bound whereas any additional nucleons will 
be loosely bound. If the energy gaps are larger than 
the strength of residual interactions such as pairing or 
quadrupole-quadrupole interactions, only weak correla- 
tions of such types are induced, hence, no corresponding 
low-lying collectivity is observed. 

Magic numbers for stable nuclei follow harmonic- 
oscillator shells up to N, Z = 20; for heavier nuclei, 
the spin-orbit splitting between high-^ orbitals produces 
N,Z = 28, 50, 82, and 126. The emergence of new 
magic numbers far from stability is thought to be re- 
lated to the effect of the tensor force Q , which produces 
(among others) attractive contributions for 0d 3 / 2 neu- 
trons in the presence of 0g?5/2 protons. For oxygen iso- 
topes the proton 0^5/2 orbit is not occupied. This con- 
tributes to the reduced binding of the v{Qd^/2) orbital 



[3| and produces a gap at N — 16, while for stable nu- 
clei the gap at N — 16 disappears since the 0d 5 / 2 proton 
shell is almost filled. The effect of the occupation num- 
ber on SP spectra is taken into account by considering 
effective SP (ESP) energies, where SP energies are mod- 
ified by the average monopole contribution of residual 
interactions. The neutron ESP spectra for neutron-rich 
oxygen isotopes derived with the universal sd-shell (USD) 
interaction (which takes into account the effects of resid- 
ual forces by virtue of a fit to experimental data) reveal 
large gaps for 22 O and 24 O at the Fermi energy. Both 
nuclei are predicted to be doubly magic with N = 14 
and 16, respectively Q; hence, the properties of 23 O are 
expected to be well reproduced by a SP description. In 
22 O, a high- lying first excited state has been observed in 
7-ray spectroscopy [H confirming the magic character of 
N = 14 for Z = 8. Systematics of neutron separation 
energies suggest a magic character of N = 16 at the neu- 
tron dripline Q , but no high-lying first excited state has 
been found for 24 O [5j. Focusing on the SP nucleus 23 O, 
we have populated the n- 22 continuum by 2p + In re- 
moval from a 26 Ne beam and deduced the decay-energy 
spectrum of excited 23 O in order to search for resonances 
corresponding to unbound states. Two-proton knockout 
reactions have proven to be a sensitive tool to study nu- 
clear structure of rare isotopes Q. These direct reac- 
tions rely on the knockout of valence protons leaving the 
remaining nucleons largely undisturbed. In the present 
paper, direct knockout of inner-shell protons is utilized 
for the first time to observe neutron-unbound states in 
neutron-rich nuclei. It is shown that these proton hole 
states couple to specific highly excited neutron states. 

The experiment was performed at the National Su- 
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FIG. 1: Experimental setup for measuring fragment-neutron 
coincidences at the NSCL. The A1900 and extended focal 
plane scintillator are upstream and not shown here. 



perconducting Cyclotron Laboratory (NSCL) at Michi- 
gan State University. A 105 pnA primary beam of 
140 MeV/A 40 Ar impinged on a 893 mg/cm 2 Be pro- 
duction target. The resulting cocktail beam was puri- 
fied with respect to the desired 26 Ne at 86 MeV/A us- 
ing an achromatic 750-mg/cm 2 -thick acrylic wedge de- 
grader and either 1% or 3% momentum-acceptance slits 
at the dispersive focal plane for different parts of the 
experiment, as well as ±10 mm slits at the extended fo- 
cal plane of the A1900 fragment separator Q. A pu- 
rity of up to 93.2% was achieved with a 26 Ne beam in- 
tensity of about 7000 pps. The contaminants (mainly 
27 Na and 29 Mg) were separated event-by-event in the off- 
line analysis by their different time-of-flight (ToF) from 
the extended focal plane to the scintillator in front of 
the 721-mg/cm 2 -thick Be reaction target (see Fig. [TJ. 
Positions and angles of incoming beam particles were 
measured by two 15 x 15 cm 2 position-sensitive parallel- 
plate avalanche counters (PPACs) with 20 pads/inch, i.e., 
FWHM « 1.3 mm. Due to the focusing of the quadrupole 
triplet, positions of impinging 26 Ne particles on the re- 
action target could be reconstructed within a radius of 
ay = 1 mm. 

Charged particles behind the reaction target were bent 
43° by the large-gap 4 Tm Sweeper Magnet Q. Two 
30 x 30 cm 2 cathode-readout drift chambers (CRDCs) 
provided position in the dispersive (10 pads/inch, a x « 
1.1 mm) and non-dispersive (drift time, a y ~ 1.3 mm) di- 
rection. The 1.87 m distance between the CRDCs trans- 
lates this into a = 1 mrad angle resolution. Energy loss 
was determined in a 65-cm-long ion chamber (IC) and a 
40 x 40 cm 2 , 4.5-mm-thin plastic scintillator whose pulse- 
height signal was corrected for position. Energy loss was 
used to separate reaction products with different Z (see 
Fig. The thin scintillator also gave ToF of reaction 
products from the reaction target. This, together with 
the total kinetic energy (TKE) measurement in a 15-cm- 
thick plastic scintillator, provided isotopic separation (see 
Fig. [2b, d). The pulse- height signal of the thick scintil- 
lator was also corrected for position; the raw ToF was 
corrected for (i) position on the thin scintillator, for (ii) 
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FIG. 2: a) Energy loss in thin scintillator vs. ion chamber 
provides Z identification, b) TKE in thick scintillator vs. 
corrected ToF for oxygen fragments provides isotope identifi- 
cation, c) Neutron ToF for decay of 23 O*, i.e., in coincidence 
with 22 O fragments, and using a thick target in singles mode 
showing prompt 7 radiation (inset). The 1.4 ns FWHM of 
the 7 peak translates into a n-energy resolution of 1.4 MeV. 
d) Data in b) projected onto the diagonal line. 



angle and position behind the Sweeper Magnet, as well 
as for (iii) the point of interaction on the target. The 
ToF corrections approximately compensate the spread in 
energy of reaction products and account for differences 
in length of the fragment tracks and of the light paths 
in the thin scintillator. Due to the finite acceptance a 
magnetic-rigidity range of reaction products is selected. 
For reaction products with equal Z, this translates into 
TKE oc 1/A and ToF oc A (see Fig. [2b). 

Neutrons were detected in the Modular Neutron Ar- 
ray (MoNA) [l(| at a distance of 8.2 m from the reaction 
target. MoNA consists of 9 x 16 stacked 2-m-long plastic 
scintillator bars which are read out on both ends by pho- 
tomultiplier tubes (PMT). The bars are mounted hor- 
izontally and perpendicular to the beam axis. Position 
along the vertical and along the beam axis is determined 
within the thickness of one bar (10 cm, a ~ 3 cm). Hor- 
izontal position and neutron ToF are determined by the 
time difference and the mean time, respectively, of the 
two PMT signals which yield resolutions of a m 5 cm and 
a ps 0.1 ns. The ToF spectrum was calibrated by shifting 
the prompt 7-ray peak from a thick-target MoNA singles 
run to 27 ns, see inset in Fig. [2b. A neutron ToF spec- 
trum (first hit after ~ 27 ns) in coincidence with 22 
fragments is given in Fig. [2b. 
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FIG. 3: Decay-energy spectra of 23 0*. The inset shows a 
close-up of the first 360 keV. Smooth curves correspond to 
simulations (see text). The solid curve corresponds to the 
sum of the dash-dotted and dashed curves. 



The decay energy of resonances is reconstructed by the 
invariant mass method. For this, the relativistic four- 
momentum vectors of the neutron and fragment are re- 
constructed at the point of breakup. For neutrons, po- 
sition and ToF resolution translate into angle and en- 
ergy resolution of a = 8 mrad and a = 1.6 MeV [lj|, 
respectively. The angle and energy of fragments in co- 
incidence with neutrons were reconstructed behind the 
reaction target based on the ion-optical properties of 
the Sweeper Magnet using a novel method which takes 
into account the position at the target in the dispersive 
direction [TjJ. The reconstructed position in the non- 
dispersive direction serves thereby as a double check on 
an event-by-event basis against the same information ob- 
tained from forward tracking from the PPACs through 
the quadrupole triplet. The angle and energy resolution 
of the fragments behind the target are a — 2.7 mrad and 
a = 0.4 MeV/^4, respectively. The average energy loss of 
the fragment through half of the reaction target is added 
to approximate the relativistic four-momentum vector of 
the fragment at the average breakup point. 

The experiment ran for two days and produced 5700 n- 
22 coincidences. The reconstructed decay-energy spec- 
trum from n- 22 coincidences is shown in Fig. [3j This 
spectrum is affected by the finite acceptance for frag- 
ments and neutrons and a bias toward the fastest neu- 
tron for events where more than one neutron hit MoNA. 
The latter is because in the analysis only the first hit 
after ~ 27 ns is taken into account. Later hits due to ei- 
ther multiple scattering of one neutron or true multiple- 
neutron events are ignored in the analysis ll| . The most 
severe of these effects, namely the acceptance cuts, have 
been simulated. For the simulation, a Glauber reac- 
tion model is used. Angle straggling of the fragments 
in the target is taken into account, as well as detec- 



tor resolutions. The n- 22 data in Fig. [3] are best de- 
scribed by a Breit-Wigner resonance at 45(2) keV decay 
energy (dashed curve) on top of a beam- velocity source of 
Maxwellian-distributed neutrons (a thermal model) with 



FIG. 4: a) Likely configuration of the 24 continuum popu- 
lated by direct 2p removal from 26 Ne. These negative-parity 
proton excitations will mix with certain negative-parity neu- 
tron excitations b), c), which then decay into the 23 O ground 
state or continuum. 



T ~ 0.7 MeV (dash-dotted curve). The ratio of the con- 
tributions is 1:2. The simulated decay-energy resolution, 
in qualitative agreement with [II], is found to scale as 
FWHM - 16V E ■ keV and AOVE ■ keV, below and above 
E = 1.5 MeV decay energy, respectively. The numerical 
factors are dominated by n-angle resolution and target 
thickness, respectively. For the present resonance, the 
width due to experimental conditions (100 keV) over- 
shadows the Wigner limit by about three orders of mag- 
nitude, hence its lifetime is not determined. The thermal 
model has been included in the fit because it is a good 
description of uncorrelated events as determined from 
event mixing, and it serves as an estimate of non-resonant 
contributions to the decay-energy spectrum (see, e.g., 
for further discussions of experimental back- 
grounds). 

When combined with the neutron separation energy 
S„ [3], our resu lt establishes an excitation energy of 
2.79(13) MeV for the first excited state of 23 0. Two low- 
lying states with spin-parities of 3/2 + and 5/2 + , based 
on a v{0d z / 2 ) 1 particle and ^(0d 5 / 2 ) _1 hole configura- 
tion, respectively, are predicted by theory [15|, with the 
5/2+ level the lower one and hence being the favored 
spin-parity assignment for the here observed resonance. 
Observation of decay to excited states in 22 O which are 
at 3.4 MeV and higher are thought to be unlikely. 

The selective population of the states can be under- 
stood from the structure of 24 O and the specific reac- 
tion mechanism. The main component (90%) of the 24 
ground-state wavefunction has the closed-shell configura- 
tion [d®, 2 , s i/2]- The 26 Ne ground-state wavefunction is 
dominated by two sd-shell protons coupled to the ground 
state of 24 Two-proton removal from 26 Ne to the 
low-lying states in 24 O can thus proceed by the removal 
of two sd-shell protons to the 24 O ground state or by one 
sd-shell and one p-shell proton to the negative-parity con- 
figuration shown in Fig. [4^,. The negative-parity eigen- 
states are linear combinations of the three configurations 
in Fig.|U The neutron-decay of these eigenstates goes by 
d 3 / 2 neutron decay to high- lying negative parity states 
(b) and by /p-shell neutron decay to the l/2 + and 5/2 + 
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one-neutron- hole states (c) in 23 O. Thus, starting with 
the main component of the 26 Ne wavefunction no overlap 
with the 3/2+ [dg / 2 , £^3/2] excited state of 23 is possible. 
The direct reaction theory for the removal of two sd- 
shell protons from the 26 Ne to the 24 O ground state gives 
a cross section of about 0.42 mb where a suppression fac- 
tor of 0.5 was used 17J. The removal of one proton from 
the sd-shell and one proton from the fully occupied p- 
shell (6 protons) corresponds to a total cross section of 
about 2.1 mb. These states will neutron decay to the 
low-lying l/2 + and 5/2 + states of 23 0. In addition, we 
might expect some direct three-nucleon knockout from 
removal of two sd-shell protons plus one sd-shell neutron, 
but again, with the dominant [rf5/ 2 ' s i/2] neu t ron struc- 
ture of 26 Ne this can only go to the l/2 + [dS, 2 , s {/2\ an d 



5/2+ [dl /21 s\ /2 ] states of 23 0. Population of the 3/2+ 
\d^ 2 ,d 3 / 2 ] state of 23 as well as the 2+ [d^ 2 ,s 1 / 2 ,d 3 / 2 ] 



state of 24 can only come from the small component 
of the 26 Ne ground state with the neutron configuration 
[g?5^ 2 , Si/2> ^3/2] < X'2 + (7r). With full sd-shell wavefunctions 
the cross section to the 2 + state of 24 is 0.03 mb (with 
the reduction factor of 0.5) compared to 0.42 mb for the 
ground state. Accordingly, the cross section for the 3/2 + 
state in 23 O should be smaller by about an order of mag- 
nitude than those for the l/2 + ground and 5/2 + excited 
states. 

Considering only the dominant components of the 
wavefunctions the spectroscopic factor for the neutron 
decay of the 5/2+ [d\f 2 ,s\^ 2 ] state in 23 to the [d\j 2 ] 
ground state in 22 O is zero. With the full USDB wave- 
functions [IH, the spectroscopic factor is 0.059. Thus, 
with a d-wave neutron single-particle width of 90(10) eV 
(corresponding to a decay energy of 45(2) keV), the total 
decay width of 5.0(6) eV is extremely small. Still, the 
calculated 7-decay lifetime is 4.5 ps corresponding to a 
partial width of 0.15 meV, hence the decay is dominated 
by neutrons. 

In conclusion, we have observed for the first time an 
excited state in 23 O. This state has been sought for un- 
successfully using 7-ray spectroscopy; due to its unbound 
nature it is revealed in the present work in an rt-fragment 
coincidence experiment. On the other hand, a predicted 
second excited state has not been found. This has been 
explained by the presence of selection effects in the pop- 
ulation of states in 23 O, which, in return, demonstrates 
that inner-shell proton knockout is a new promising spec- 
troscopic tool to explore excited states far from stability. 
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The nominal energy resolution of a — 0.7 MeV is de- 
graded by a systematic ToF walk with deposited energy. 
This walk does neither affect the position resolution along 
a bar, nor can it be corrected due to inefficiencies of the 
charge-signal collection. Furthermore, the effect on the 
decay energy resolution is negligible as it is determined 
mostly by the neutron angle resolution and the target 
thickness. 



